Spatial transcriptomics reveals altered lipid metabolism and inflammation-related gene expression of sebaceous glands in psoriasis and atopic dermatitis.

Sebaceous glands drive acne, however, their role in other inflammatory skin diseases remains unclear. To shed light on their potential contribution to disease development, we investigated the spatial transcriptome of sebaceous glands in psoriasis and atopic dermatitis patients across lesional and non-lesional human skin samples. Both atopic dermatitis and psoriasis sebaceous glands expressed genes encoding key proteins for lipid metabolism and transport such as ALOX15B, APOC1, FABP7, FADS1/2, FASN, PPARG, and RARRES1. Also, inflammation-related SAA1 was identified as a common spatially variable gene. In atopic dermatitis, genes mainly related to lipid metabolism (e.g. ACAD8, FADS6, or EBP) as well as disease-specific genes, i.e., Th2 inflammation-related lipid-regulating HSD3B1 were differentially expressed. On the contrary, in psoriasis, more inflammation-related spatially variable genes (e.g. SERPINF1, FKBP5, IFIT1/3, DDX58) were identified. Other psoriasis-specific enriched pathways included lipid metabolism (e.g. ACOT4, S1PR3), keratinization (e.g. LCE5A, KRT5/7/16), neutrophil degranulation, and antimicrobial peptides (e.g. LTF, DEFB4A, S100A7-9). In conclusion, our results show that sebaceous glands contribute to skin homeostasis with a cell type-specific lipid metabolism, which is influenced by the inflammatory microenvironment. These findings further support that sebaceous glands are not bystanders in inflammatory skin diseases, but can actively and differentially modulate inflammation in a disease-specific manner.

Generative artificial intelligence empowers digital twins in drug discovery and clinical trials.

INTRODUCTION: The concept of Digital Twins (DTs) translated to drug development and clinical trials describes virtual representations of systems of various complexities, ranging from individual cells to entire humans, and enables in silico simulations and experiments. DTs increase the efficiency of drug discovery and development by digitalizing processes associated with high economic, ethical, or social burden. The impact is multifaceted: DT models sharpen disease understanding, support biomarker discovery and accelerate drug development, thus advancing precision medicine. One way to realize DTs is by generative artificial intelligence (AI), a cutting-edge technology that enables the creation of novel, realistic and complex data with desired properties. AREAS COVERED: The authors provide a brief introduction to generative AI and describe how it facilitates the modeling of DTs. In addition, they compare existing implementations of generative AI for DTs in drug discovery and clinical trials. Finally, they discuss technical and regulatory challenges that should be addressed before DTs can transform drug discovery and clinical trials. EXPERT OPINION: The current state of DTs in drug discovery and clinical trials does not exploit the entire power of generative AI yet and is limited to simulation of a small number of characteristics. Nonetheless, generative AI has the potential to transform the field by leveraging recent developments in deep learning and customizing models for the needs of scientists, physicians and patients.

The Oncology Biomarker Discovery framework reveals cetuximab and bevacizumab response patterns in metastatic colorectal cancer.

Precision medicine has revolutionised cancer treatments; however, actionable biomarkers remain scarce. To address this, we develop the Oncology Biomarker Discovery (OncoBird) framework for analysing the molecular and biomarker landscape of randomised controlled clinical trials. OncoBird identifies biomarkers based on single genes or mutually exclusive genetic alterations in isolation or in the context of tumour subtypes, and finally, assesses predictive components by their treatment interactions. Here, we utilise the open-label, randomised phase III trial (FIRE-3, AIO KRK-0306) in metastatic colorectal carcinoma patients, who received either cetuximab or bevacizumab in combination with 5-fluorouracil, folinic acid and irinotecan (FOLFIRI). We systematically identify five biomarkers with predictive components, e.g., patients with tumours that carry chr20q amplifications or lack mutually exclusive ERK signalling mutations benefited from cetuximab compared to bevacizumab. In summary, OncoBird characterises the molecular landscape and outlines actionable biomarkers, which generalises to any molecularly characterised randomised controlled trial.

The impact of the cardiovascular component and somatic mutations on ageing.

Mechanistic insight into ageing may empower prolonging the lifespan of humans; however, a complete understanding of this process is still lacking despite a plethora of ageing theories. In order to address this, we investigated the association of lifespan with eight phenotypic traits, that is, litter size, body mass, female and male sexual maturity, somatic mutation, heart, respiratory, and metabolic rate. In support of the somatic mutation theory, we analysed 15 mammalian species and their whole-genome sequencing deriving somatic mutation rate, which displayed the strongest negative correlation with lifespan. All remaining phenotypic traits showed almost equivalent strong associations across this mammalian cohort, however, resting heart rate explained additional variance in lifespan. Integrating somatic mutation and resting heart rate boosted the prediction of lifespan, thus highlighting that resting heart rate may either directly influence lifespan, or represents an epiphenomenon for additional lower-level mechanisms, for example, metabolic rate, that are associated with lifespan.

Oncogene-induced MALT1 protease activity drives posttranscriptional gene expression in malignant lymphomas.

Constitutive mucosa-associated lymphoid tissue lymphoma translocation protein 1 (MALT1) activity drives survival of malignant lymphomas addicted to chronic B-cell receptor signaling, oncogenic CARD11, or the API2-MALT1 (also BIRC3::MALT1) fusion oncoprotein. Although MALT1 scaffolding induces NF-κB-dependent survival signaling, MALT1 protease function is thought to augment NF-κB activation by cleaving signaling mediators and transcriptional regulators in B-cell lymphomas. However, the pathological role of MALT1 protease function in lymphomagenesis is not well understood. Here, we show that TRAF6 controls MALT1-dependent activation of NF-κB transcriptional responses but is dispensable for MALT1 protease activation driven by oncogenic CARD11. To uncouple enzymatic and nonenzymatic functions of MALT1, we analyzed TRAF6-dependent and -independent as well as MALT1 protease-dependent gene expression profiles downstream of oncogenic CARD11 and API2-MALT1. The data suggest that by cleaving and inactivating the RNA binding proteins Regnase-1 and Roquin-1/2, MALT1 protease induces posttranscriptional upregulation of many genes including NFKBIZ/IκBζ, NFKBID/IκBNS, and ZC3H12A/Regnase-1 in activated B-cell-like diffuse large B-cell lymphoma (ABC DLBCL). We demonstrate that oncogene-driven MALT1 activity in ABC DLBCL cells regulates NFKBIZ and NFKBID induction on an mRNA level via releasing a brake imposed by Regnase-1 and Roquin-1/2. Furthermore, MALT1 protease drives posttranscriptional gene induction in the context of the API2-MALT1 fusion created by the recurrent t(11;18)(q21;q21) translocation in MALT lymphoma. Thus, MALT1 paracaspase acts as a bifurcation point for enhancing transcriptional and posttranscriptional gene expression in malignant lymphomas. Moreover, the identification of MALT1 protease-selective target genes provides specific biomarkers for the clinical evaluation of MALT1 inhibitors.

Numerical and Machine Learning Analysis of the Parameters Affecting the Regionally Delivered Nasal Dose of Nano- and Micro-Sized Aerosolized Drugs.

The nasal epithelium is an important target for drug delivery to the nose and secondary organs such as the brain via the olfactory bulb. For both topical and brain delivery, the targeting of specific nasal regions such as the olfactory epithelium (brain) is essential, yet challenging. In this study, a numerical model was developed to predict the regional dose as mass per surface area (for an inhaled mass of 2.5 mg), which is the biologically most relevant dose metric for drug delivery in the respiratory system. The role of aerosol diameter (particle diameter: 1 nm to 30 µm) and inhalation flow rate (4, 15 and 30 L/min) in optimal drug delivery to the vestibule, nasal valve, olfactory and nasopharynx is assessed. To obtain the highest doses in the olfactory region, we suggest aerosols with a diameter of 20 µm and a medium inlet air flow rate of 15 L/min. High deposition on the olfactory epithelium was also observed for nanoparticles below 1 nm, as was high residence time (slow flow rate of 4 L/min), but the very low mass of 1 nm nanoparticles is prohibitive for most therapeutic applications. Moreover, high flow rates (30 L/min) and larger micro-aerosols lead to highest doses in the vestibule and nasal valve regions. On the other hand, the highest drug doses in the nasopharynx are observed for nano-aerosol (1 nm) and fine microparticles (1-20 µm) with a relatively weak dependence on flow rate. Furthermore, using the 45 different inhalation scenarios generated by numerical models, different machine learning models with five-fold cross-validation are trained to predict the delivered dose and avoid partial differential equation solvers for future predictions. Random forest and gradient boosting models resulted in R2 scores of 0.89 and 0.96, respectively. The aerosol diameter and region of interest are the most important features affecting delivered dose, with an approximate importance of 42% and 47%, respectively.

Machine Learning and Pharmacometrics for Prediction of Pharmacokinetic Data: Differences, Similarities and Challenges Illustrated with Rifampicin.

Pharmacometrics (PM) and machine learning (ML) are both valuable for drug development to characterize pharmacokinetics (PK) and pharmacodynamics (PD). Pharmacokinetic/pharmacodynamic (PKPD) analysis using PM provides mechanistic insight into biological processes but is time- and labor-intensive. In contrast, ML models are much quicker trained, but offer less mechanistic insights. The opportunity of using ML predictions of drug PK as input for a PKPD model could strongly accelerate analysis efforts. Here exemplified by rifampicin, a widely used antibiotic, we explore the ability of different ML algorithms to predict drug PK. Based on simulated data, we trained linear regressions (LASSO), Gradient Boosting Machines, XGBoost and Random Forest to predict the plasma concentration-time series and rifampicin area under the concentration-versus-time curve from 0-24 h (AUC0-24h) after repeated dosing. XGBoost performed best for prediction of the entire PK series (R2: 0.84, root mean square error (RMSE): 6.9 mg/L, mean absolute error (MAE): 4.0 mg/L) for the scenario with the largest data size. For AUC0-24h prediction, LASSO showed the highest performance (R2: 0.97, RMSE: 29.1 h·mg/L, MAE: 18.8 h·mg/L). Increasing the number of plasma concentrations per patient (0, 2 or 6 concentrations per occasion) improved model performance. For example, for AUC0-24h prediction using LASSO, the R2 was 0.41, 0.69 and 0.97 when using predictors only (no plasma concentrations), 2 or 6 plasma concentrations per occasion as input, respectively. Run times for the ML models ranged from 1.0 s to 8 min, while the run time for the PM model was more than 3 h. Furthermore, building a PM model is more time- and labor-intensive compared with ML. ML predictions of drug PK could thus be used as input into a PKPD model, enabling time-efficient analysis.

Infection Control Measures and Prevalence of SARS-CoV-2 IgG among 4,554 University Hospital Employees, Munich, Germany.

Hospital staff are at high risk for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection during the coronavirus disease (COVID-19) pandemic. This cross-sectional study aimed to determine the prevalence of SARS-CoV-2 infection in hospital staff at the University Hospital rechts der Isar in Munich, Germany, and identify modulating factors. Overall seroprevalence of SARS-CoV-2-IgG in 4,554 participants was 2.4%. Staff engaged in direct patient care, including those working in COVID-19 units, had a similar probability of being seropositive as non-patient-facing staff. Increased probability of infection was observed in staff reporting interactions with SARS-CoV-2‒infected coworkers or private contacts or exposure to COVID-19 patients without appropriate personal protective equipment. Analysis of spatiotemporal trajectories identified that distinct hotspots for SARS-CoV-2‒positive staff and patients only partially overlap. Patient-facing work in a healthcare facility during the SARS-CoV-2 pandemic might be safe as long as adequate personal protective equipment is used and infection prevention practices are followed inside and outside the hospital.

Multi-omic landscaping of human midbrains identifies disease-relevant molecular targets and pathways in advanced-stage Parkinson's disease.

BACKGROUND: Parkinson's disease (PD) is the second most common neurodegenerative disorder whose prevalence is rapidly increasing worldwide. The molecular mechanisms underpinning the pathophysiology of sporadic PD remain incompletely understood. Therefore, causative therapies are still elusive. To obtain a more integrative view of disease-mediated alterations, we investigated the molecular landscape of PD in human post-mortem midbrains, a region that is highly affected during the disease process. METHODS: Tissue from 19 PD patients and 12 controls were obtained from the Parkinson's UK Brain Bank and subjected to multi-omic analyses: small and total RNA sequencing was performed on an Illumina's HiSeq4000, while proteomics experiments were performed in a hybrid triple quadrupole-time of flight mass spectrometer (TripleTOF5600+) following quantitative sequential window acquisition of all theoretical mass spectra. Differential expression analyses were performed with customized frameworks based on DESeq2 (for RNA sequencing) and with Perseus v.1.5.6.0 (for proteomics). Custom pipelines in R were used for integrative studies. RESULTS: Our analyses revealed multiple deregulated molecular targets linked to known disease mechanisms in PD as well as to novel processes. We have identified and experimentally validated (quantitative real-time polymerase chain reaction/western blotting) several PD-deregulated molecular candidates, including miR-539-3p, miR-376a-5p, miR-218-5p and miR-369-3p, the valid miRNA-mRNA interacting pairs miR-218-5p/RAB6C and miR-369-3p/GTF2H3, as well as multiple proteins, such as CHI3L1, HSPA1B, FNIP2 and TH. Vertical integration of multi-omic analyses allowed validating disease-mediated alterations across different molecular layers. Next to the identification of individual molecular targets in all explored omics layers, functional annotation of differentially expressed molecules showed an enrichment of pathways related to neuroinflammation, mitochondrial dysfunction and defects in synaptic function. CONCLUSIONS: This comprehensive assessment of PD-affected and control human midbrains revealed multiple molecular targets and networks that are relevant to the disease mechanism of advanced PD. The integrative analyses of multiple omics layers underscore the importance of neuroinflammation, immune response activation, mitochondrial and synaptic dysfunction as putative therapeutic targets for advanced PD.

Inferred Ancestral Origin of Cancer Cell Lines Associates with Differential Drug Response.

Disparities between risk, treatment outcomes and survival rates in cancer patients across the world may be attributed to socioeconomic factors. In addition, the role of ancestry is frequently discussed. In preclinical studies, high-throughput drug screens in cancer cell lines have empowered the identification of clinically relevant molecular biomarkers of drug sensitivity; however, the genetic ancestry from tissue donors has been largely neglected in this setting. In order to address this, here, we show that the inferred ancestry of cancer cell lines is conserved and may impact drug response in patients as a predictive covariate in high-throughput drug screens. We found that there are differential drug responses between European and East Asian ancestries, especially when treated with PI3K/mTOR inhibitors. Our finding emphasizes a new angle in precision medicine, as cancer intervention strategies should consider the germline landscape, thereby reducing the failure rate of clinical trials.

Artificial intelligence in early drug discovery enabling precision medicine.

Introduction: Precision medicine is the concept of treating diseases based on environmental factors, lifestyles, and molecular profiles of patients. This approach has been found to increase success rates of clinical trials and accelerate drug approvals. However, current precision medicine applications in early drug discovery use only a handful of molecular biomarkers to make decisions, whilst clinics gear up to capture the full molecular landscape of patients in the near future. This deep multi-omics characterization demands new analysis strategies to identify appropriate treatment regimens, which we envision will be pioneered by artificial intelligence.Areas covered: In this review, the authors discuss the current state of drug discovery in precision medicine and present our vision of how artificial intelligence will impact biomarker discovery and drug design.Expert opinion: Precision medicine is expected to revolutionize modern medicine; however, its traditional form is focusing on a few biomarkers, thus not equipped to leverage the full power of molecular landscapes. For learning how the development of drugs can be tailored to the heterogeneity of patients across their molecular profiles, artificial intelligence algorithms are the next frontier in precision medicine and will enable a fully personalized approach in drug design, and thus ultimately impacting clinical practice.

A statistical framework for assessing pharmacological responses and biomarkers using uncertainty estimates.

High-throughput testing of drugs across molecular-characterised cell lines can identify candidate treatments and discover biomarkers. However, the cells' response to a drug is typically quantified by a summary statistic from a best-fit dose-response curve, whilst neglecting the uncertainty of the curve fit and the potential variability in the raw readouts. Here, we model the experimental variance using Gaussian Processes, and subsequently, leverage uncertainty estimates to identify associated biomarkers with a new Bayesian framework. Applied to in vitro screening data on 265 compounds across 1074 cancer cell lines, our models identified 24 clinically established drug-response biomarkers, and provided evidence for six novel biomarkers by accounting for association with low uncertainty. We validated our uncertainty estimates with an additional drug screen of 26 drugs, 10 cell lines with 8 to 9 replicates. Our method is applicable to any dose-response data without replicates, and improves biomarker discovery for precision medicine.

Identification of Intrinsic Drug Resistance and Its Biomarkers in High-Throughput Pharmacogenomic and CRISPR Screens.

High-throughput drug screens in cancer cell lines test compounds at low concentrations, thereby enabling the identification of drug-sensitivity biomarkers, while resistance biomarkers remain underexplored. Dissecting meaningful drug responses at high concentrations is challenging due to cytotoxicity, i.e., off-target effects, thus limiting resistance biomarker discovery to frequently mutated cancer genes. To address this, we interrogate subpopulations carrying sensitivity biomarkers and consecutively investigate unexpectedly resistant (UNRES) cell lines for unique genetic alterations that may drive resistance. By analyzing the GDSC and CTRP datasets, we find 53 and 35 UNRES cases, respectively. For 24 and 28 of them, we highlight putative resistance biomarkers. We find clinically relevant cases such as EGFRT790M mutation in NCI-H1975 or PTEN loss in NCI-H1650 cells, in lung adenocarcinoma treated with EGFR inhibitors. Interrogating the underpinnings of drug resistance with publicly available CRISPR phenotypic assays assists in prioritizing resistance drivers, offering hypotheses for drug combinations.

Concise Polygenic Models for Cancer-Specific Identification of Drug-Sensitive Tumors from Their Multi-Omics Profiles.

In silico models to predict which tumors will respond to a given drug are necessary for Precision Oncology. However, predictive models are only available for a handful of cases (each case being a given drug acting on tumors of a specific cancer type). A way to generate predictive models for the remaining cases is with suitable machine learning algorithms that are yet to be applied to existing in vitro pharmacogenomics datasets. Here, we apply XGBoost integrated with a stringent feature selection approach, which is an algorithm that is advantageous for these high-dimensional problems. Thus, we identified and validated 118 predictive models for 62 drugs across five cancer types by exploiting four molecular profiles (sequence mutations, copy-number alterations, gene expression, and DNA methylation). Predictive models were found in each cancer type and with every molecular profile. On average, no omics profile or cancer type obtained models with higher predictive accuracy than the rest. However, within a given cancer type, some molecular profiles were overrepresented among predictive models. For instance, CNA profiles were predictive in breast invasive carcinoma (BRCA) cell lines, but not in small cell lung cancer (SCLC) cell lines where gene expression (GEX) and DNA methylation profiles were the most predictive. Lastly, we identified the best XGBoost model per cancer type and analyzed their selected features. For each model, some of the genes in the selected list had already been found to be individually linked to the response to that drug, providing additional evidence of the usefulness of these models and the merits of the feature selection scheme.

Stratification and prediction of drug synergy based on target functional similarity.

Drug combinations can expand therapeutic options and address cancer's resistance. However, the combinatorial space is enormous precluding its systematic exploration. Therefore, synergy prediction strategies are essential. We here present an approach to prioritise drug combinations in high-throughput screens and to stratify synergistic responses. At the core of our approach is the observation that the likelihood of synergy increases when targeting proteins with either strong functional similarity or dissimilarity. We estimate the similarity applying a multitask machine learning approach to basal gene expression and response to single drugs. We tested 7 protein target pairs (representing 29 combinations) and predicted their synergies in 33 breast cancer cell lines. In addition, we experimentally validated predicted synergy of the BRAF/insulin receptor combination (Dabrafenib/BMS-754807) in 48 colorectal cancer cell lines. We anticipate that our approaches can be used for prioritization of drug combinations in large scale screenings, and to maximize the efficacy of drugs already known to induce synergy, ultimately enabling patient stratification.

Quantitative Proteome Landscape of the NCI-60 Cancer Cell Lines.

Here we describe a proteomic data resource for the NCI-60 cell lines generated by pressure cycling technology and SWATH mass spectrometry. We developed the DIA-expert software to curate and visualize the SWATH data, leading to reproducible detection of over 3,100 SwissProt proteotypic proteins and systematic quantification of pathway activities. Stoichiometric relationships of interacting proteins for DNA replication, repair, the chromatin remodeling NuRD complex, ß-catenin, RNA metabolism, and prefoldins are more evident than that at the mRNA level. The data are available in CellMiner (discover.nci.nih.gov/cellminercdb and discover.nci.nih.gov/cellminer), allowing casual users to test hypotheses and perform integrative, cross-database analyses of multi-omic drug response correlations for over 20,000 drugs. We demonstrate the value of proteome data in predicting drug response for over 240 clinically relevant chemotherapeutic and targeted therapies. In summary, we present a novel proteome resource for the NCI-60, together with relevant software tools, and demonstrate the benefit of proteome analyses.

Defining subpopulations of differential drug response to reveal novel target populations.

Personalised medicine has predominantly focused on genetically altered cancer genes that stratify drug responses, but there is a need to objectively evaluate differential pharmacology patterns at a subpopulation level. Here, we introduce an approach based on unsupervised machine learning to compare the pharmacological response relationships between 327 pairs of cancer therapies. This approach integrated multiple measures of response to identify subpopulations that react differently to inhibitors of the same or different targets to understand mechanisms of resistance and pathway cross-talk. MEK, BRAF, and PI3K inhibitors were shown to be effective as combination therapies for particular BRAF mutant subpopulations. A systematic analysis of preclinical data for a failed phase III trial of selumetinib combined with docetaxel in lung cancer suggests potential indications in pancreatic and colorectal cancers with KRAS mutation. This data-informed study exemplifies a method for stratified medicine to identify novel cancer subpopulations, their genetic biomarkers, and effective drug combinations.